Interference of Instanton Trajectories in Quantum Tunneling for Small Particles of 

Real Antiferromagnets 
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For a two-sublattice antiferromagnet the Lagrangian is constructed taking into account Berry 
phase whose form is matched with the quantum-mechanical Heisenberg Hamiltonian. Tunnel effects 
are analyzed taking into account the crystallographic symmetry and possible types of Dzyaloshinskii- 
Moriya interaction. It is shown that, when the real magnetic symmetry and the Dzyaloshinskii- 
Moriya interaction are taken into consideration, the effects of a destructive instanton interference 
and the suppression of macroscopic quantum tunneling can play an essential role. It also may lead 
to a periodic dependence of the ground-state level splitting on the Dzyaloshinskii-Moriya interaction 
constant; the magnitude of this splitting is calculated. 
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I. INTRODUCTION 

During the last decade, macroscopic quantum tunnel- 
ing in macroscopic (or, to be more precise, mesoscopic) 
magnetic systems has become an object of intense exper- 
imental and theoretical investigations.^ In the physics of 
magnetism such systems include small magnetic parti- 
cles, magnetic clusters, and high-spin molecules. Special 
attention is paid to the coherent macroscopic quantum 
tunneling (CMQT) between physically different, but en- 
ergetically equivalent, states in systems with discrete de- 
generacy of the ground state. A typical CMQT effect in 
such systems is tunneling between two equivalent classi- 
cal states corresponding to two minima of the anisotropy 
energy.^ 

The CMQT effects can be observed experimentally 
from the resonant absorption of electromagnetic waves 
at tunnel-splitted energy levels. The interest in these ef- 
fects are associated with the two following factors. First, 
mesoscopic objects exhibiting quantum-mechanical prop- 
erties are interesting as potential elements for quantum 
computers. Second, fine and elegant effects of interfer- 
ence of instanton trajectories emerge in these problems. 
For ferromagnetic particles these effects suppress tun- 
neling for half-integral values of the total spin of the 
system'^"'* and lead to oscillatory dependences of the tun- 
nel splitting of energy levels on extrinsic parameters.^ In 
addition, in contrast to the effects of quantum escape 
from a metastable to stable state the manifestations of 
the CMQT effects are not masked by thermal fluctua- 
tions. 

Initially, the CMQT investigations were carried out for 
small particles of a ferromagnet®'^ under the assumption 
that all spins in the particle are parallel to one another 
(the high-spin model). The effects of destructive inter- 
ference of instanton trajectories and interference sup- 
pression of tunneling were predicted precisely for such 
systems. ■^'"''^ It turned out later that antiferromagnets 
form a more convenient class for experimental investiga- 
tions of CMQT. According to calculations of Refs.^'^°, 
the level spitting in antiferromagnets is stronger than in 



ferromagnets, and the effects can be observed at high 
temperatures. It is not surprising that the CMQT effects 
were observed for the first time in ferritin particles with 
an antiferromagnetic structure.^ No interference effects 
are observed in pure antiferromagnets (i.e., in the case of 
complete compensation of the spins of sublattices), but 
such effects may appear in the applied magnetic field. ^ 
It will be shown below that even at zero field the inter- 
ference effects can also appear when the real magnetic 
symmetry of the crystal is taken into account, in partic- 
ular, in the presence of the Dzyaloshinskii-Moriya (DM) 
interaction. 

A semiclassical description of magnetic systems is 
based on the formalism of coherent spin states. In order 
to construct the effective field Lagrangian both for ferro- 
magnets and antiferromagnets, we will proceed from the 
expression for the Euclidean Lagrangian of an individual 
spin, which has the form^^ 



(1) 



Here, s is the spin associated with each magnetic mo- 
ment, (j)k and 9k are the polar coordinates of the kth 
magnetic moment, and W{4'k,0k) is the classical mag- 
netic energy of the magnet; the overdot indicates the 
differentiation with respect to the imaginary time t = it. 
The first term determines the magnetization dynamics 
(its variation leads to the well-known Landau-Lifshitz 
equations in the angular parametrization) and also de- 
termines a so-called Berry phase. ^^'^^ This quantity is 
associated with the total time derivative which is not 
manifested itself in the equations of motion, but it is re- 
sponsible for the interference of instanton trajectories. 

For a macroscopic description it is natural to use one 
or several field variables (order parameters) instead of 
the set of microscopic variables. The determination of 
the number of order parameters and their transformation 
properties in magnetic systems is a nontrivial problem. 
In the approach based on coherent spin states the order 
parameter for ferromagnets is the magnetization vector 
of a constant length, which can be parametrized by the 
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angular variables 6 and (j). In this case the Berry phase is 
just a resultant change in the angle along the instanton 
trajectory.^'^ The behavior of antiferromagnetic systems 
can be correctly described using a three-component vec- 
tor of a fixed length, viz., the antiferromagnetism vector 
/, see Refs.^'^""'^^. The total spin in this case is a slave vari- 
able, and it is determined by the vector I and its time 
derivative dl/dt. Dynamic equations for the antiferro- 
magnetism vector I can be either constructed from the 
symmetry considerations^'* or derived from the Landau- 
Lifshitz equations for the sublattice magnetizations.*^'*^ 
In both these approaches the same classical equations of 
motion for the unit vector I are obtained, which are usu- 
ally referred to as the equations of the u-model. The 
application of such equations considerably simplifies the 
analysis of both linear and nonlinear dynamic effects in 
an antiferromagnet, see Refs.*^ *^. However, the advan- 
tage of these equations for describing macroscopic quan- 
tum effects is not so obvious. The Lagrangian obtained 
from the classical Landau-Lifshitz equations or from sym- 
metry considerations cannot be used directly for describ- 
ing the MQT effect taking into account the interference 
of instanton trajectories. It is probably for this reason 
that Golyshev and Popkov^ used in their analysis of the 
CMQT effects a system of two equations for the sublat- 
tice magnetizations, whose analysis is much more com- 
plicated. 

As a matter of fact, it is impossible in principle to re- 
construct the Lagrangian of a dynamic system from the 
classical equations of motion. The Lagrangians describ- 
ing the same classical equations of motion for the system 
can differ in the term which is the total derivative with 
respect to time. This term does not effect on the clas- 
sical dynamics of the system, but alters the magnitude 
of the Euclidean action on trajectories. For this rea- 
son the corresponding terms with total derivatives were 
lost in the early publications.^'^ A consistent qiiantum- 
mechanical expression for the spin Lagrangian taking into 
account the correct equation for the total derivative can 
be derived using the formalism of coherent states and the 
analysis of the evolution operator; this expression coin- 
cides with Eq. (1) ibove. Topological terms of the 
form of total derivatives in the effective Lagrangian for 
the vector I are significant for the quantum theory of ID 
antiferromagnets.** However, it is impossible in principle 
to derive their expressions only from the classical equa- 
tions of the (T-model for the vector I. 

In the simplest version of the u-model the derivatives 
of I with respect to time appear in the Lagrangian in 
the trivial form {dl/dr)"^, see Refs.*^~*^. In this case 
the equations of the ct- model are Lorcntz-invariant, and 
the description of the dynamics of nonlinear magneti- 
zation waves (kink-type solitons in antiferromagnets) is 
considerably simplified.*'''*^ The interference effects in 
the MQT are obviously absent. However, the situation 
changes drastically for more realistic models. 

First, for many antiferromagnetic crystals there ex- 
ist terms reflecting interactions of the DM type, which 



are linear in I and in magnetization. It was shown in 
Refs.^°'^* that these interactions are responsible for the 
terms in the effective Lagrangian which are linear in 
dl/dr; this considerably modifies the kink dynamics in 
comparison with the simplest Lorentz-invariant model. 
Obviously, such interactions can also lead to the emer- 
gence of total derivatives (topological phases). The pres- 
ence of a magnetic field may also lead to similar effects; 
this was noted in the analysis of the nonlinear dynamics 
of antiferromagnets^^ as well as for the MQT effect, see 
Ref.^^ and recent publications^'^^"^'^. 

In the present paper we will construct the Lagrangian 
of the (T-model on the base of the Eq. (1) taking into 
account all possible sources of the terms with the to- 
tal derivative, which may lead to nontrivial interference 
effects. This Lagrangian will be used to study the in- 
terference of instanton trajectories for the real models of 
antiferromagnetic particles of various symmetries and to 
determine the contribution of these effects to the tunnel- 
ing probability. 

II. LAGRANGIAN OF THE cr-MODEL FOR 
REAL ANTIFERROMAGNETS 

Let us consider a system with localized spins, in which 
nearest neighbors are coupled through the antiferromag- 
netic interaction. We assume that the lattice has such 
a structure that the sites with spins can be divided into 
two groups so that the spins appearing in pairs of near- 
est neighbors belong to different groups and there are 
no frustrations in the lattice. For ideal antiferromagnets 
these two groups correspond to the two magnetic sublat- 
tices. In this case we assume that the spins correspond- 
ing to each group have parallel orientations and form 
the total spins Si and ^2 of the sublattices. In the ex- 
change approximation for such antiferromagnets, vectors 
Si and <S'2 are antiparallel. The total spin S = Si + S2 
in the ground state can differ from zero in view of a dif- 
ferent number of sites in the sublattices (decompensa- 
tion), I Si I ^ 1, and also in the presence of an external 
magnetic field and/or the DM interaction, when the an- 
tiparallelism of the spins is violated (i.e., \Si + S2\ ^ 
even for = IS2I). We will consider only completely 
compensated antiferromagnets with |S'i| = |S'2|, since 
the specific effects associated with spin decompensation 
( I Si I ^ 1^2 1, but the difference |S'i-S'2| < |Si,2|) reduce 
the interference effects to those which are well known for 
ferromagnets.^* 

Our goal is to construct the Lagrangian describing the 
dynamics of the vector I in the presence of the DM in- 
teraction and the magnetic field. Since the tensor of 
exchange interaction constants Jy may have the anti- 
symmetric component in the nearest neighbor approxi- 
mation, the Hamiltonian of such a system has the form 

ne = J J2 ^-'^f^+ J2 d-{S„xS0)-gtiBj2^-^»- 

<afS> <oi(3> a 

(2) 
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Here, the first term describes the isotropic exchange in- 
teraction, the summation in this term is extended to the 
pairs of nearest neighbors, and Sa is the spin at the ath 
site. The antisymmetric component of the tensor of the 
exchange constants is a microscopic source of the DM 
interaction,^^ c;orrcsponding to the dual vector d. The 
last term describes the interaction of spins with the ex- 
ternal magnetic field. 

Let us consider the exchange approximation in which 
the deviation from the conventional Heisenberg model 
with an isotropic exchange interaction JSaSjs is small; 
i.e., d, gjisH <C J. In this case we can introduce the total 
spins 5^1=^ Sa^ and 3^ = ^ Sa^ of the sublattices and 
assume that the vectors Si and S2 have a fixed length. 
It is convenient to put Si = Nsa\ and ^2 = Nscr2, 
where s is the spin at a site and A'' is the number of sites 
in each sublattice. We will parametrize the unit vectors 
CTi and (72 by the polar coordinates {6\,(j)\) and {02, 4>2), 
respectively. In this case the classical magnetic energy 
of the antiferromagnet, whose exchange component cor- 
responds to the Hamiltonian (2), can be written in the 
form 

W{ai, (72) = Js^zNcri(T2 + s^zNd ■ (cri x a2) 

+ w{cri,(T2) - gnBaNH{ai+ CT2) ■ (3) 

Here, N is the number of spins in sublattices, z is the 
coordination number for a lattice site, and w{cr\,cr2) is 
the anisotropy energy. 

Thus, we arrive at the description of the energy of an 
antiferromagnet in terms of two unit vectors. Their dy- 
namics can be described by a Lagrangian which can be 
written in the dynamic variables <t\ and cr2 taking into 
account relation (3) in the form 

C = -ihSiAi{(Ti)di-ihS2A2{cr2)(T2-W{(Ti,CT2) . (4) 

Here, we have chosen a more general form of the kinetic 
terms as compared to relation (1). These terms can be 
presented through the vector potential of the field of a 
magnetic monopole: 



^1,2(0-) 



cr X ni,2 
a{a + crni,2) 



(5) 



where ni^2 are the quantization axes of coherent states 
for each sublattice. This potential has a singularity for 
<7n = —a, i.e., on a certain half-line in the cr-space. 
Usually, the "north pole" gauge with n = is used, in 
which the quantity A{cr)& assumes to have the famil- 
iar form (1). The potentials Ai^2 of the monopole field 
permit gauge transformations (such as a change in the 
position of spin quantization axes and, hence, singulari- 
ties), which do not change the equations of motion, but 
make a contribution to the Lagrangian in the form of the 
total derivative of the function of spins (Xi and a2 with 
respect to r, which can be significant for the description 
of interference effects. The kinetic terms for each sub- 
lattice can be written in individual gauges (in particular. 



with different directions of the quantization axes ni and 
n2)- 

We present the unit vectors cri and CT2 in terms of 
vectors I = (<7i — cr2)/2 and m = {cri +0-2) /2 which are 
connected through the relations 



I' 



= 1 



and 



(6) 



In the (T-model approximation, when the magnetic mo- 
ment is small \m\ <^ \l\, we can use simple transforma- 
tions in order to present the Lagrangian (4) as a power 
expansion in m. Confining the expansion to the terms 
linear in m, we can write the kinetic terms in the form 



ihAi{<Ti)di - ihA2{(T2)d2 
-ihi[Ai{l) - A2{-1)] - ihrrii 



dA,{l) ^ ■dA2{-l) 



+ 1- 



dli dli 
- ihm[Ai{l) + A2{-1)] . (7) 

Here, the form of the vector-potentials Ai and A2 has 
not been specified yet. In particular, the quantization 
axes n\ ans Ti2 have not been chosen. It is natural to 
choose the quantization axes n\ and n2 so that A\{1) = 
A2{—1), which is possible for ni = — n2- In this case 
the singular term with dl/dr vanishes, and the dynamic 
terms expansion starts with the term linear in m, which 
can be written in the form 



—ihm ■ {F X I) — ih—[m ■ Ai{l)] , 
dr 



where 



Fi = e 



ijk 



fdA, dAk 



dli 



(8) 



(9) 



Thus, most of the arbitrariness in the choice of the gauge 
field A, which takes place for ferromagnets. does not ex- 
ist for antiferromagnets. The gauge-invariant quantity 
Fj, which has the meaning of a formal magnetic field 
associated with potential A, is the magnetic monopole 
field F = In the transition from Eq. (7) to the 

expression (8) the initial gauge arbitrariness turned out 
to be localized in the term with the total derivative 
d[mA{l)]/dT. Concerning this quantity, its contribution 
to the Euclidean action is obviously equal to zero in the 
case when an instanton trajectory misses the singular 
point of A{1). This condition can be easily satisfied if 
we choose the direction n = ni = —n2 along the hard 
magnetization axis of the antiferromagnet. In this case 
the phase for a closed path on the unit sphere = 1, 
which is formed by instanton trajectories, is independent 
of the position of the quantization axis n. 

Taking into account the condition ml = 0, we elimi- 
nate from the expression (4) the slave variable m: 



m 



2Jsz 



-fiH""^ - liH''^ ■ I)) - il X / 



(10) 



where 7 = giig/h is the gyromagnetic ratio and H'^^ is 
the effective field which is the sum of the external field 
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H and the DM field ff^M ^y^c approximation chosen 
above, in whieh the DM interaction can be presented in 
the purely antisymmetric form d{cri x 0-2) 3c d{l x m), 
the DM field can be written as H^^ = zs{d x l)/{gi2B)- 
The expression for m is also valid for more general 
forms of the DM interaction, which cannot be reduced 
to a bilinear form in <Ti^2- In particular, we will consider 
more general forms of the DM interaction of the type 
Dik{l)milk which arc observed for many crystals and are 
significant for the MQT effects. In this case the effective 
field in the expression for m assumes the form Hf^ = 
Dik{l)lk, and 



Hf = hI°^ + Dik{l)lk ■ 



(11) 



In this section we will not specify the form Dik{l). The 
approximation \m\ <C \l\ used in the derivation of the a- 
model is satisfied for max(_ff, H^^) <C H^x, where i?ex = 
Jsz/ jiB is the exchange field. Substituting m into the 
Lagrangian (4), we obtain the effective Lagrangian for 
the vector I in the form 



C = 



2Jz 



1 • 



+ i^H^*" -{Ixl) 



Wail) 



l{H'''^'-l)-H 



DM 



(12) 



Here, Wa{l) has the meaning of the effective anisotropy 
energy in which the additional term is (H^^ ' 0^ ~ 
(^fj'DM^2 ^;a,ken into account along with the initial energy 
w{l) = w{(Ti,cr2) introduced above for cri = —cr^ = I- 
The quantity Wa{l) is obviously the real anisotropy en- 
ergy determined from static measurements in weak fields, 
and there is no point in separating these contributions. 
We must simply use the expression for Wa{l) which is 
determined by the crystal symmetry of the magnet. The 
specific form of the anisotropy energy for various anti- 
ferromagnets is given in the Table I. The terms in the 
braces describe the variation of the static energy of the 
antiferromagnet due to the external magnetic field. The 
first term, which is quadratic in the components of H, 
is quadratic in I and can also be presented as the field 
induced renormalization of the anisotropy energy. The 
second term, which is bilinear in the components of the 
external magnetic field H and the DM field H^^'^, con- 
tains odd powers of the components of I and describes the 
energy of the weak ferromagnetic moment induced by the 
DM interaction. (In particular, this term can be reduced 
to H-{dx I) for a purely antisymmetric DM interaction.) 
This term can completely remove the degeneracy of the 
classical ground state of the system, and the analysis of 
the MQT effects becomes meaningless. For this reason, 
it makes sense to take into account the external magnetic 
field and the DM interaction simultaneously only for cer- 
tain selected orientations of the external field, when this 
term vanishes for vector I directed along the easy magne- 
tization axis of the antiferromagnet. Some of these ori- 



entations of the field for orthorhombic antiferromagnets 
were considered in Ref.^. 

Thus, we arrive at the following conclusions. The La- 
grangian for the vector I differs from the Lagrangian for 
the (T-modcl of ideal antiferromagnets^^ in the presence 
of a number of additional terms which play different roles 
in the description of CMQT. In contrast to the case of fer- 
romagnets or antiferromagnets with different spins of the 
sublattices, the term with the total derivative can easily 
be eliminated. It is important to note that the external 
field and some forms of the DM interaction leads to the 
emergence of gyroscopic terms linear in dl/dr. The emer- 
gence of these terms indicates the lowering of the actual 
dynamic symmetry of antiferromagnets in the presence 
of a magnetic field and/or the DM interaction. 

The structure of the Lagrangian is such that the con- 
tribution of the DM interaction to the gyroscopic term 
can be taken into account by adding the DM field H^^, 
which is a function of Z, to the external magnetic field H. 
Gyroscopic terms can make significant contributions to 
the probabilities of tunneling processes both by effecting 
on the structure of instanton solutions and by creating 
destructive interference of instanton trajectories. It will 
be proved below that, in contrast to the case of a ferro- 
magnets or antiferromagnets with different spins of the 
sublattices, this interference is not of topological origin, 
but can also be given below. The examples of "pure" 
antiferromagnets in which tunneling can be completely 
suppressed due to the interference of instanton trajecto- 
ries will be given below. 

Deriving the Lagrangian (12), we neglected the possi- 
bility of inhomogeneous tunneling and, hence, the depen- 
dence of S and I on spatial coordinates was omitted from 
the very outset. The inclusion of such a dependence leads 
to the substitution N ^ J dV/a^, where a is the lattice 
constant, and to the emergence of an additional term pro- 
portional to Ja^iyi)^ in the Euclidean action. A com- 
parison of the inhomogeneity energy with the anisotropy 
energy leads to an estimate of the spatia l inhomogeneity 
size on the order of Aq = Hex/ Han, where H^x and 
Han are the exchange field and the anisotropy field and 
Aq is the domain wall thickness. If the size of a particle is 
larger than Aq, i.e., iV > TV, ~ (Ao/a)^ ~ {H^x/Hanf^, 
we can assume a more advantageous inhomogeneous tun- 
neling scenario, in which the level splitting weakly de- 
pends on N (or is even independent of it) for N > Nc 
Although this question has not been discussed in the lit- 
erature and its analysis is beyond the scope of the present 
publication, we will briefly consider it. 

The value of Nc is too large for the tunneling effects to 
be observable for N > Nc. As a matter of fact, N in the 
tunneling exponent is multiplied by the susceptibility of 
the system, i.e., appears in the combination N Han /Hex, 
see Ref."'^''. The presence of this small parameter actually 
makes it possible to observe tunneling at ferritin parti- 
cles with N ~ 3.5-10^, see Refs.^'^°. However, for typical 
values of Han/ Hex ~ 10~^ - 10~^, the tunneling expo- 
nent NcHan/Hex ^ [Hex/HanY^"^ » 1 is too large and 
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the observation of the transition to the inhomogeneous 
tunneling mode becomes problematic. 



III. SYMMETRY OF THE INSTANTON 
SOLUTIONS AND INTERFERENCE OF 
CONTRIBUTIONS FROM INSTANTON 
TRAJECTORIES 

In accordance with the general rules of the semiclassi- 
cal approximation formulated in the instanton language, 
the amplitude of transition from one state to another is 
described in the so-called instanton-gas approximation^^. 
The level splitting for a system with two equivalent min- 
ima can be presented in the form 



A oc D\/K , 
where the quantity D is defined as 



D = {det'n) 



1/2 



V 2Trh ) 



exp 



ReX \ 



(13) 



(14) 



and T is the one-instanton action; K is a, combinatorial 

factor emerging due to nonuniqueness of the tunnel path 
connecting two equivalent minima; and det 'fl is the fluc- 
tuation determinant disregarding the zeroth mode, which 
is determined by small deviations from an instanton tra- 
jectory, see Ref.''^ for details. In order to analyze the ef- 
fects of tunneling between degenerate states correspond- 
ing to the ground states of the system and to determine 
the value of splitting, we must find the one-instanton 
trajectories connecting these states, calculate the value 
of the Euclidean action I on these trajectories, and find 
the determinant of the operator for the second variation 
of action. The contribution to the splitting comes only 
from equivalent trajectories corresponding to the mini- 
mum value of the real component of X. The combinato- 
rial factor depending on the phase difference in the tra- 
jectories will be calculated below using Eq. (18). In this 
section, we concentrate our attention on an analysis of 
the main contribution which comes only from X and will 
not calculate the preexponential factor. Let us see how 
these calculations can be carried out in actual practice. 

For a concrete analysis, it is convenient to write the 
Lagrangian in the form 



C = 



2Jz 



^fdiy ., „ dl 

-(-] +^iu:,xl).- + fMl) 



(15) 



where u>h = jH'^^, 7 is the gyromagnetic ratio, and H'^^ 
is the effective field. The dimensionless function Wa(0 is 
proportional to the anisotropy energy, and the value of 
luq coincides with the frequency of a homogeneous anti- 
ferromagnetic resonance in the uniaxial anisotropy field. 
We parametrize the vector I by the angular variables 



h = sin 6 cos c 



I2 = sin 6 sin ( 



I3 = cos 9 . 



We are dealing with an easy-axis anisotropy. Conse- 
quently, the ground state is doubly degenerated and has 
two values of I corresponding to it: I = 63 and I = —63, 
and the unit vector 63 being parallel to the easy axis. 
Let us consider the tunneling between these two states. 
Function Wa{i) for a magnet with the anisotropy axis C„ 
can be written in the form 



(17) 



where the first term corresponds to easy-axis anisotropy 
and Wa{0, 0) <C 1 defines anisotropy in the basal plane. 

For antiferromagnets with an easy axis of symmetry 
C„, there exist n instanton trajectories and n antiin- 
stanton trajectories, and the combinatorial factor has the 
form 



K 



k,k'=0 



(16) 



^fe,S' = Tlm/ dT£{lj), 
Jk\Jk' 

(18) 

i.e., ^/ is the phase difference between the kth instan- 
ton and the k'th antiinstanton. The integral defining 
j^/ is taken over a closed path formed by the trajec- 
tories of the fcth instanton and the k'th antiinstanton. 
In the Lorentz-invariant a-model, the Lagrangian is real 
and all ^f. j., are equal to zero; the combinatorial factor 

K = 'n? \s trivial and equal to n. Consequently, ^/K = n; 
i.e., the total transition amplitude and level splitting for 
nn pairs is just the contribution from one instanton mul- 
tiplied by the number of paths. It will be shown below, 
however, that for j^, 7^ 0, the level splitting A may 
contain an oscillatory dependence on the product of the 
small parameter |u>_ff| and the large quantity N, and, 
hence, requires a more detailed analysis. The nature of 
its oscillations can be established from symmetry consid- 
erations, and the specific form of the function K of the 
parameters of the problem can be determined even with- 
out solving the corresponding the Euler-Lagrange equa- 
tions. 



A. Lorentz-invariant cr-model 

It is convenient to consider first the tunneling in the 
simplest Lorentz-invariant cr-model which corresponds to 
the Lagrangian (15) with ljh = 0. As a matter of fact, 
for some models of an antiferromagnet with the DM in- 
teraction, the results turn out to be the same as in the 
absence of this interaction, see below. If H'^^ = 0, the 
analysis of the problem does not present any difficulty. 
Indeed, for any form of the anisotropy energy in a uniax- 
ial antiferromagnet with the principal axis C^2 , ^4, (in 
the subsequent analysis, we will consider only the type 
of symmetry that can exist in the crystal lattice), the 
instanton solution corresponds to the function 9 — 9(t) 
with the boundary conditions ^ 0, tt for r ^ ±00 and 
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<j) = <j)Q = const, where is defined by the relation 

dwa{9,(t>) 



(19) 



Let us assume that the ground states ±63 are on the 
principal axis C„. In this case the value of w'a{9,4') is 
proportional to sin ncj) and there exist 2n solutions to this 
equation: 



ttA; 



fc = 0, 1,.. .,2n - 1 , 



(20) 



from which n solutions 



fe,mm 



correspond to the minima 

of Wa{9,(f)), while the remaining n solutions (pfl^^^^ cor- 
respond to the maxima of this function for all (9 7^ 0, tt. 



Instantons with 



/c,min 



correspond to the lowest value of 



the Euclidean action, and we will consider below only 
these n solutions. Function 9{t) can be determined from 
the second-order equation for which the first integral is 
known to be 



^oK(0,0i°V^«a(O,</.i°))]. (21) 



Henceforth, we assume that Wa{0,(j)) = Wa{7T,(j)) = 
and that the value of = corresponds to the minimum 
of the function Wa{9, (j)). With such a choice of the axes, z 
is always an easy magnetization axis and x is an medium 
magnetization axis. The Euclidean action on trajectories 
is real- valued for all values of 6 and is defined as 



J = 



2Jz 



d9^wa{9,<j>) . 



(22) 



This approximate expression is written in the main 
approximation in small anisotropy in the basal plane /3 ^ 
1 , where [3 is the characteristic anisotropy constant in the 
basal plane, i.e., the maximum value of Wa ■ Thus, the 
contribution in the given case comes from n instanton 
trajectories on which vector I is real and rotates in one of 

the n planes defined by the condition (j) = 4>^k^^i-a = ■ 
The imaginary component of T in the Lorentz-invariant 
model is absent, and the combinatorial factor K in the 
expression (18) is equal to n^. 



B. Role of if^" 

The inclusion of the terms with H'^^, which destroy the 
Lorentz invariance, brings about two types of difficulties. 
First, for H'^^ ^ 0, the solution </) = const is generally 
inapplicable, and the instanton structure is determined 
by the general system of two second-order equations 



)sm 



sin 6' cos 6' + Wq^^ + iu>H^^{9, 

Uu 



sm f cos ( 



■LUr, 



iujH0T{9, 



0, (23a) 

-0, 
(23b) 



whose solutions are generally not real-valued. Here, the 
terms with F are determined by the variation of the term 
with • [{dl/dr) x I] in the Lagrangian (15), and the 
form of the function r{9,(f>) generated by the DM inter- 
action for various types of magnetic symmetry is given 
in the column 5 of the Table I. Second, the imaginary 
component of the Euclidean action X, which comes from 
the term proportional to (jJh, may appear even for tra- 
jectories with a real I. Let us consider the cases when 
these situations are realized. 

If F(^, (/)) vanishes at the same values of 0^,°^ as for 
dwa{9, (f>)/d(t>, the second equation (23b) is satisfied iden- 
tically for the plane trajectories </> = 0, while the first 
equation in the system (23a) can be reduced to Eq. (21) 
considered above in the Lorentz-invariant cr-model. Con- 
sequently, in this case T{9, (p) does not effect on the form 
of the function 9 = 6{t) in the instanton solution, but 
changes the imaginary component of the Euclidean ac- 
tion. This effect will be considered in more details in the 
section IV. 

If, however, T{9,(j)l^^) ^ 0, the instanton does not cor- 
respond to a plane solution </> = const any longer, and we 
must seek the general solution of the system (23) in the 
form 9 = 9{t), (j) = (j){T). In this case the functions 9{t) 
and (P{t) may in general turn out to be complex-valued. 
There are no general analytical methods for constructing 
such separatrix solutions; an instanton solution of the 
system of equations (23) can be written exactly only for 
some cases (see Ref.^^ and the section IV in the present 
paper). 

It will be shown below that the effect of the term in the 
Lagrangian on the imaginary component of the Euclidean 
action I may lead to nontrivial consequences even for 
antiferromagnets with T{e,(f)) ^ 0, but F(6l,(/)^°^) = 0, 
and there exists a real- valued instanton solution 6 = 6{t), 

(p = or in the case when the value of H'^^/H^x is 

negligibly small or its inclusion changes 9 = 9{t) and the 
real component of I insignificantly. 

In order to explain this fact, we consider the case when 
the value of T /Hf.x <C 1 is so small that instanton tra- 
jectories can be regarded as planar, 9 = 9{t), (j) = const. 
The presence of the term linear in dl /dr leads to the con- 
tribution to the imaginary component of the Euclidean 
action I, which is proportional to the number of spins 
in the particle. The imaginary component of the Eu- 
clidean action I is of order of ImX/fi, cx Nd/J: i.e., it 
is proportional to the product of a small and a large 
parameter. Consequently, the effects of destructive in- 
terference can be significant. It is well known that the 
interference effects for orthorhombic ferromagnets may 
suppress tunneling completely."'^"'* In contrast to the case 
of antiferromagnets the term with dm/dr for ferromag- 
nets does not contain a small factor H^^/Hex, but it is 
inessential since the value of Im I/ft ~ nNs ^ 1 for ferro- 
magnets, while tunneling is completely suppressed when 
Im J/?l ~ TT. This condition can easily be satisfied for a 
large TV. In particular, for the antiferromagnetic particle 
of the ferritin with N ~ 3500 the tunneling probability 
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TABLE I: Anisotropy in the basal plane and the Dzyaloshinskii-Moriya interaction constant for systems with various types of 
the magnetic symmetry. 



1 


2 


3 




5 


6 


7 


8 


n 


Wa 


Axes 


JJMl 


1 (6*, (p) 




, 0) 








^2 ^3: 


TTlx^y H~ TTlylx 


3 sin^ 61 sin 20 


* 


3 sin^ sin 20 


4 


2 


P2 sm f) Sin <p 


Z,2 Zf3j 




6 sin^ cos ^ sin 


* 


6 sin cos sin 


4 






9(-)9(-)9(+) 

Zg; Zy 


rUxlz + TWj/a! 


6 sin^ 6 cos cos cj) 




6 sin ^ cos 6 cos 


4 








i(m+4-m-Z?.) 


5 sin® ^ sin 40 


* 


5 sin* sin 40 


16 


4 


Bi sin'' 6 sin^ 2S 


4(-)2(+)2(-) 


777 ... / — Tn^.I^. 

' ' "X ''X ' ' ''y^y 


3 sin^ S cos 2(A 




3 sin"^ cos 20 


16 






4(-)2(-)2(+) 

^z ^xy 




3sin^6lsin20 


* 


3 sin^ sin 20 


8 + 8cos(sdAr/J) 






6i^^2i-)2t7) 




7 sin'^ sin 60 


* 


7 sin® 61 sin 60 


36 


6 


fSe sin« 9 sin^ 30 


6<-'2<+'2l7) 




5 sin"' d cos 6 sin 30 


* 


5 sin^ 6 cos sin 30 


36 








hriAi^'- +/') 


5 sill* cos cos 'id 




5 siii'^ cos cos 3r') 


36 



"For high-order axes, the following notation is introduced: m± = 
mx ± iniy and i± = ± iZy. 

'Asterisks mark systems for which an exact solution correspond- 
ing to the munimum of the real part of the action exist. 



in the magnetic field with interference is an osciUating 
function of the field, and the suppression of tunneling 
can be observed in fields H < 100 Oe, see Refs.^'^®'^^, 
which are much weaker than the characteristic value of 
the DM field H^^ = 10^ - 10^ Oe. 

On the other hand, the contribution to the real com- 
ponent of the Euclidean action does not contain the large 
parameter A''. This contribution can be appreciable, see 
the next section, but in this case the product of other 
parameters, namely, the small quantity d/J <C 1 and 
the large quantity » 1, is significant. Thus, the 
terms with dl/dr may lead to two types of effects: (i) 
the emergence of nonplanar instanton trajectories and 
complex values of components of I on these trajectories; 
(ii) the interference of instantons even in the case of plane 
trajectories with the real Euclidean action I. 

The first effect only takes place when the term r{6, 4>) 
in Eq. (23) differs from zero. Such terms arc always im- 
portant for the description of the domain wall dynamics 
in antiferromagnets: they may reduce the limit veloc- 
ity of the domain wall motion to a considerable extent 
and may also lead to an abrupt change in the wall struc- 
ture upon a continuous variation of its velocity.^"'^^ The 
subsequent analysis of concrete instanton solutions will 
show that the role of such terms in the description of 
the instanton structure and tunneling is not so impor- 
tant as in the description of the domain wall dynamics. 
On the other hand, if the function T{9, 0) differs from 

zero, but the function = for the given solution 

V{9,(j)'^^) = 0, the domain wall dynamics is trivial and 
can be described by Lorentz-invariant expressions. In 



this case, the instanton structure 9 — 9{t) is the same as 
in the Lorentz-invariant theory. However, the situation 
with instantons is different: not all features can be de- 
scribed by the function 6{t) and the real component of 
the Euclidean action T only. It will be shown bellow that 
the main contribution from the term H'^^ ■ [{dl/dT) x I] 
is associated precisely with interference processes and is 
manifested most clearly exactly when the instanton tra- 
jectory is planar; i.e., r{9, 0*-°'') = 0. 

In the case of real-valued trajectories, it is convenient 
to use the following approach for calculating the imagi- 
nary component of the Euclidean action. '^^ We introduce 
the vector r = rl which is not subjected to the condition 
r"^ — 1 and present the term with the first derivative in 
the expression (15) in the form 

-i^A— , with A = 5 . (24) 

This expression has the same structure as the term in 
the nonrclativistic Lagrangian describing the interaction 
of a classical charged particle moving in a 3D space with 
coordinate r and velocity v = dr/dr with a formal mag- 
netic field B = V X A (differentiation is carried out in 
the r space). It is well known that the magnetic field 
appears in the Lagrangian of a charged particle through 
the vector potential A at point r, which is defined only 
to within a certain gauge, while the field B is gauge- 
invariant. 

Simple but cumbersome calculations proved that, for 
any ferromagnet, this formal magnetic field B may be 
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radial and can be presented in the form 

B = ^B{6,<I>), (25) 

where 

B{e, <!>) = 2(ff^ff/) - ^ + . (26) 

Oli oik 

In the absence of the DM interaction the value of B(0, </>) 
is determined only by the external field H^^\ B{6,4>) = 
2{H''°h). At zero field the value of 5(6', (j)) is determined 
by the DM field Hf^^ = Dij{l)lj and can be presented 
in terms of the tensor Dij as 

B{0, (j)) = SDijlilj — Dii + Dij^kkljlk — Dij^ilj . (27) 

Here, the comma in the subscript in D indicates the dif- 
ferentiation of the tensor Dij with respect to the corre- 
sponding component of /, and the summation is extended 
over double indices. The values of B{9, 0) for various 
types of DM interaction and of the configurations of axes 
are given in column 7 of the Table I. 

Phases ^k,k' defined in Eq. (18) can be presented 
in terms of the integrals jA.dr taken over instanton- 
antiinstanton pairs forming closed loops. Using the 
Stokes theorem, we can present the phase diff'erence ^k,k' 
as the magnetic flux of the vector B through a part of 
the unit sphere bounded by such a loop. Obviously, indi- 
vidual phases are determined by the vector potential 
i.e., depend on the gauge, but the phase differences are 
gauge-invariant for all pairs of trajectories. 

It is important that the structure of B{9, (f) for all 
possible types of DM interaction is such that the total 
flux of field B through the unit sphere, 

/■TT /'27r 

$t„t=/ smOdO I d(i>B{6,(j)) (28) 
Jo Jo 

is equal to zero and the value of cos is indepen- 
dent of the total time derivative in the Lagrangian (gauge 
invariance).^^ 

This approach enables us to calculate speciflcally the 
phase difi^erence of integral trajectories and the combina- 
torial factor K in Eq. (13). We begin with the simplest 
case of an orthorhombic antiferromagnet for which there 
are only two pairs of instanton trajectories. It can be 
easily verified that the DM interaction of any type (see 
the Table I) makes a zero contribution to the phase dif- 
ference for two diametrically opposite trajectories.^* For 
this reason the required phase factor can be determined 
only by the external field and can be written in the form 

cos$ = cosl cosa 1 , (29) 

where a is the angle between the plane containing the 
instanton trajectories and the external field iJ^"^. This 
result was obtained by Chiolero and Loss^'^ in particular 



cases a = and a = tt/2. Thus, for instanton trajectories 
lying in the same plane all possible types of DM interac- 
tion given in the Table I do not effect on the tunneling. 
This results does not contradict the analysis carried out 
by Golyshev and Popkov^, who studied tunneling in small 
completely compensated particles of an orthorhombic an- 
tiferromagnet with the orthoferrite structure and discov- 
ered that no interference effects exist in a zero magnetic 
field. The our approach enabled us to obtain this result 
without solution of the Euler-Lagrange equations. Thus, 
we have proved that the conclusion that no interference 
takes place for diametrically opposite trajectories can be 
extended to more general cases of the DM interaction. 
It is important that this conclusion is drawn without re- 
sorting to any approximation, which is inevitable in the 
solution of a complex system of equations describing the 
instanton structure. 

Thus, in the case of orthorhombic antiferromagnets 
with two instanton trajectories none of the types of DM 
interaction presented in the Table I leads to destructive 
interference. However, this result is different for uniax- 
ial antiferromagnets with an easy magnetization axis C„, 
n > 2. In this case there exist n pairs of instanton trajec- 
tories. Obviously, here we also have pairs of trajectories 
lying in the same plane, for which k' ~ ^ ^^"^ inter- 
ference is trivial, but it can be manifested itself for pairs 
of trajectories with 0^°^ — 0^?^ ^ tt. 

It will be shown below that the value of combinatorial 
factor can be reduced considerably from its maximum 
value to zero; i.e., both partial and complete suppres- 
sion of interference are possible. 

C. Tetragonal antiferromagnets 

Let us demonstrate it using specific examples of par- 
ticles with a tetragonal easy magnetization axis and bi- 
nary axes in the perpendicular plane (crystallographic 
class 4222,22,2,), when the minimum of the real compo- 
nent of the Euclidean action corresponds to four instan- 
ton and four antiinstanton trajectories. In order to de- 
scribe tetragonal antiferromagnets, we choose the po- 
lar axis 63 along the tetragonal easy magnetization axis 
4^. Anisotropy in the basal plane is determined by the 
fourth-order invariant, see the Table I. We assume that 
(34 > 0; i.e., instanton trajectories correspond to the ro- 
tation of I in the equivalent planes zx and zy. Depending 
on the magnetic parity of the principal axis (according 
to Turov^^) and of the binary axes 2^;, 2y, or 2,xyt 
basically different types of behavior can be observed. We 
will consider them separately. 

1. System ii+hi'hi'y'' 

With such a structure of axes in an antiferromagnet 

only the antisymmetric invariant d{mxly — mylx) is usu- 
ally considered, which can be obtained from the antisym- 
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metric component of the tensor of exchange constants. 
The value of d is of order of \/JJj, see Ref.^^. This invari- 
ant determines the weak isotropic ferromagnetic moment 
when / is oriented in the basal plane. However, it is of no 
interest to us since it can be reduced to the total deriva- 
tive in the Lagrangian and, hence, gives T{0, (f)) = in 
the equations of motion and B{6,(j)) = in the imag- 
inary part of the Euclidean action. In addition, there 
exist a number of invariants of the rclativistic origin, 
which give a nonzero value of r(0, </>). The simplest of 
these invariants has the form 2(Z^ — ly){rnxly + nriylx) 
which coincides (except for the total derivative) with 
the invariant (to+^+ — m-l^_)/(2i) presented in the ta- 
ble. It can easily be seen, however, that in this case 
T{9,(j)"^^) = 0, and the instanton solution has the form 

6 = 0{t), (p = = 7rA;/2 for integer k. The value of 
B{9, (j)) is such that 



B oc sin 9 sin 40 



(30) 



and, hence, all phases <l>j, j., are equal to zero. An anal- 
ysis of other invariants, e.g., iKirixly — iriylx), leads to 
the same result (namely, the DM interaction does not ef- 
fect on tunneling in any way). This result is apparently 
independent of the model and is determined only by the 
type of magnetic symmetry. The model independence for 
dynamic effects in domain walls has the same origin, i.e., 
the DM interaction; it was demonstrated in Refs.^°'^^. 
Thus, the case of an even principal axis 41^'' may serve 
as an example that nonzero terms which are linear in 
dl/dr and cannot be reduced to the total derivative are 
not appeared in the separatrix solution and do not effect 
on tunneling in any way. 

A different situation takes place for antiferromagnets 
with an odd principal axis 4i . In this case two cases arc 
possible: when the intermediate anisotropy axes through 
which tunneling takes place are odd and when these axes 
are even. 



It is an oscillating function of the DM interaction con- 
stant d, and it takes the values in the range from to 16. 
For realistic values of N of order of 10^ - 10^ the period 
is not large; the value of AH^^/H^^ ~ IQ-^ 10"! 
for characteristic values of ~ 10'* Oe and H^x — 

10^ Oe. Since the value of the DM field is very sensi- 
tive to extrinsic parameters (e.g., the pressure or the ad- 
dition of a small amount of impurities to the crystal), 
these oscillations can be observed and monitored. An 
additional opportunity for observing interference effects 
appears when the magnetic field is taken into considera- 
tion. 



System 4 hi+hij 



Such a symmetry group is typical for the extensively 
studied weak antiferromagnet MnF2, see Ref.^*. In this 
case T{9,(j)) cx sin^ 6i cos 2(/) and r(6',(/>) ^ for all val- 



ues of 



corresponding to the minimum of the 



anisotropy in the basal plane and describing instanton 
trajectories for d — Q. For c? 7^ instanton solutions can- 
not be written in the simple form 9 = 9{t), </> ~ ■nk/2, 
fc = 0, 1, 2, 3. On the other hand, if we assume that the 
value of d is very small, we can easily find, applying the 
approximation of planar rotation and using the formula 
B{9,4>) oc sin^ cos 2(/), that the difference in the imagi- 
nary components of T for pairs of trajectories lying in the 
same plane as well as for adjacent instanton trajectories 
is equal to zero, and no interference effects take place. We 
will consider the solution for this case in the next section 
and prove that these simple regularities are preserved in 
a more rigorous analysis also when we do not require that 
<j) = 7rfc/2. We will also consider a general mechanisms 
of the timncling in the case when T{e, <i)^^.J ^ and the 
instanton solution is not real-valued. 



D. Hexagonal antiferromagnets 



2. System 4~'^2i~hi+'^ 



In this case T{e,(t>) = and ^ = 0^°' = 7rA;/2. The 
presence of the DM interaction does not effect on the 
instanton trajectories with the minimum action, which 
correspond to </> = = 9{t). However, in contrast 
to the system with 41"'''' , the contribution of the DM in- 
teraction is significant for calculating the combinatorial 
factor K in Eq. (13). It can be seen from the explicit ex- 
pression B{9,4>) oc sin^0sin2(/) that the phase difference 
for adjacent trajectories (with (/> = cj)^^^ cj) = (/>^'2i) differs 
from zero. 

Thus, the phase factor for the tunneling probability is 
given by 



K = 16 sin- 



(31) 



Let us briefly consider the case of a hexagonal princi- 
pal axis. Here, we again have three cases presented in 

the Table I. For a system with an even principal axis 
Gz^^'^x 's^^p there exists the invariant nixly — ruylx and 
the DM interaction, which cannot be reduced to a to- 
tal derivative, appears only in the fifth order in I. The 
analysis is similar to that for the system 4i^''2i ''2iy\ 
In this case also nonzero terms which arc linear in dl/dr 
and cannot be reduced to a total derivative do not effect 
on tunneling in any way. It can be verified that such 
a behavior is the same as for antiferromagnets with the 
even principal axis ni'*'^ . 

For systems with an odd principal axis, the DM inter- 
action is a cubic function of I, but it makes a zero con- 
tribution to the imaginary component of the Euclidean 
action under the assumption that instanton trajectories 
are real and planar. In the system 6i ^2^ ■'2^^g minimal 
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instanton trajectories have an imaginary component, but 
this only changes the real part of the Euclidean action. 
The analysis of this system is similar to that which will be 
carried out in the section IV for the system 4^ ■'2i^-'2ij,''. 
Consequently, the combinatorial factor K for all three 
cases has the maximum value equal to = 36. 



functions /(r), Oir) and docs not contradict this system 
and leads to the following system of two equations for 
functions / and 6: 

— O — f^smO cos 6* + Wg sin 9 cos 

-Z{-lfujDjsm^0cosh.2f = Q, (33a) 



IV. INSTANTON SOLUTION FOR 
ANTIFERROMAGNETS WITH THE 
SYMMETRY 4"*2i+^2i^* 



It was noted above that in the case of antiferromag- 
nets with the symmetry 4^ ^2\^^2xy there is no exact 
solution of the type (p = 7rfc/2, 6 = 9{t) for trajectories 
with the rotation of I in the vicinity of the easy plane 

c± 7rfc/2, and we have to analyze the complete system of 
two equations (23a) and (23b). The situation is compli- 
cated further since these equations have complex-valued 
coefficients and, in general, their complex solutions of 
the type = 9i{t) + 16*2 (r), (j) = (/'i(t) + i(j)2{T) must be 
considered. As a result, system (23) is equivalent to a dy- 
namic system with four degrees of freedom, and it is not 
integrable. No general method exists for an analysis of 
such systems. However, a comprehensive analysis can be 
carried out in the given case as well as for magnets with 
other types of symmetry, which are listed in the Table I. 

In the case of antiferromagnets with the odd tetragonal 
axis, the equations for the angular variables 6 and <p have 
the form 

-6 + ^'^ sin 9cos9 + col sin 9 cos 9(1 + 04 sin^ 9 sin^ 2(f)) 
+ 3iujD<j) sin^ 6* cos 2^ = , (32a) 



— (f) sin^ 9 — 2^9 sin 9 cos 9 + PaJq sin^ 9 sin 2<j) cos 2(j) 
-SiuJoO sin^ 9cos2(l) = 0. (32b) 

The quantity ujq defines the height of the potential bar- 
rier through which tunneling occurs, and /34 is the di- 
mensionless anisotropy constant in the basal plane. We 
assume that /34 > 0, which corresponds to instantons 
of the Lorentz- invariant a- model (ojd = 0) that passing 
through the even axes a; or y; /34 <i; 1 corresponds to 
the easy-axis limit, and ujd is proportional to the DM 
interaction constant, ujd = 7|-H'^'^| = zd/h. 

It can be easily seen that this system has the exact 
solution (j) = n{2k + l)/4, 6 = 9{t), which has been 
considered in the previous section. It determines tun- 
neling for f34 < 0, but in the case of /34 > we are 
interested in now, it corresponds to the rotation of I in 
the hard planes, does not ensure that the real part of 
the Euclidean action has a minimum, and makes a zero 
contribution to the tunneling amplitude in the instanton 
approximation. The exact solution <p = Trk/2 does not 
exist in this case. It can be seen, however, that the sub- 
stitution (j) = Trk/2 + if{T) and 9 = 9{t) with the real 



/ sin^ d + 2f9 sin 9 cos 9- /?4Wo sin^ sinh 2/ cosh 2/ 

+ 3(-l)^a;i5^sin^(9cosh2/ = . (33b) 

In addition, such a substitution renders the Lagrangian 
real-valued: 



2Jz 



- + 



{-ifujDf cosh 2/ cos 9{2 + sin^ 9) 

- i^oPi . 4,n ■ 1,2 Of 

— — sm 9 smh 2/ 



+ -sin^ 



(34) 



It can be seen that the imaginary component / — Im (j) 
of the instanton solution effects only on the real part of 
the Euclidean action. The system (33) is equivalent to 
a mechanical system with two degrees of freedom. Only 
one first integral is known for it: 



S = 



2Jz 



P 



— sm"' ( 



+ 



'^0/^4 . 4 a • u2of 

— — sm 9 smh 2/ 



(35) 



Note that f = for the separatrix solutions we are in- 
terested in. For this reason this system cannot be an- 
alyzed exactly. However, an approximate solution can 
be constructed in the physically interesting case, when 
(jJd <^ Wq, Pi ^ 1 and for any relation between cjd 
and u>o04- In order to find such a solution, we note 
that Eq. (33a) from the system (33) is transformed into 
{9)"^ = uJq siv? 9 in the zeroth approximation in the small 
parameters wd/wq and ^4. In this case the constant so- 
lution f = fo = const satisfies Eq. (33b) and gives 



sinh 2/0 = 



/340J0 



(36) 



It should be noted that the value of /o is determined by 
the ratio of two small parameters and can be appreciable. 
Using this fact, we can write a refined equation for 9{t) 



■ ojq sm 9 cos ( 



l-^sin^e 



(37) 



The approximate solution constructed by us is valid if 
^ ~ Wo sin^, i.e. 



n LUn 

Pi 



(38) 
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This condition may also hold in the case when the value 
of sinh2/o = ujd/Pa^o is of order of unity, but it is still 
violated for /34 ^ 0. In this case the situation is similar to 
that observed for the domain wall dynamics:^^ the limit 
velocity of a domain wall in a tetragonal antiferromagnet 
with an odd axis with </> ^ const tends to zero as /J4 — * 0, 
and no dynamic solution exists for (3i = Q. 

The value of the Euclidean action for the obtained so- 
lution is real and is defined by the expression 



X=^(l 
Jz 



3/34^2 y • 



(39) 



In the range of applicability of the constructed solution, 
i.e., when ojd and /?4, are small and when inequality (38) 
is satisfied, the inclusion of the DM interaction leads only 
to a small correction to the real part of the Euclidean 
action, the imaginary part being identically equal to zero. 

For models of antifcrroniagncts with binary and hexag- 
onal symmetry axes such an approximate solution can- 
not be constructed, but the analysis of these models is 
even simpler than in the case of antiferromagnets with a 
tetragonal symmetry axis. In these cases we can also ver- 
ify that, if the exact solution (j) = 4>^k^ = 27rfc/n, 9 = 9{t) 
does not exist, the solution has the same form as before: 



\- l.f{T) 

n 



e = e{T) 



(40) 



and the term in the Lagrangian which is linear in dl / dr 
contributes only to the real part of Euclidean action J. 
It can be proved, however, that function /(r) is anti- 
symmetric and proportional to the parameter uj d / {I3ojq) , 
which is always small (in contrast to the case of a tetrago- 
nal antiferromagnet, in which there appears the parame- 
ter ijju/ {I3aU](i) whose value may be large). Consequently, 
we can seek the function / <C 1 using the same pertur- 
bation theory as for the domain wall dynamics in such 
antiferromagnets. 2^ As a result, we obtain the following 
expression for the real part of the Euclidean action: 



Jz 



1 



(41) 



where the numerical factor ^ is of order of unity. Thus, 
the correction to the result typical for the Lorentz- 
invariant model is always small. Note that no interference 
effects take place in this case. We arrive at the conclusion 
that the DM interaction effects on the tunnel probability 
in hexagonal and orthorhombic antiferromagnets. 



V. CONCLUDING REMARKS 

The analysis of antiferromagnctic particles with a 
tetragonal symmetry axis shows that three possible types 
of the effect of the DM interaction on tunneling processes 
exist. The investigation of the remaining cases impor- 
tant for the analysis of crystalline antiferromagnets (or- 
thorhombic or uniaxial with a hexagonal symmetry axis. 



sec the Table I) proved that these types include all pos- 
sible cases for antiferromagnctic systems with a doubly 
degenerate ground state. In fact, all cases can be reduced 
to the following three types of behavior. 



1. The principal axis is even, e.g., 41'*'^ or 6i^\ 



In 



this case vector I is real on all instanton trajecto- 
ries and these trajectories are planar (j) = 27rfc/n, 
= 0(t). The real part of the Euclidean action 
is independent of the DM interaction constant, the 
imaginary part is equal to zero, and destructive in- 
terference effects are absent. In this case tunneling 
can be described without taking into account the 
DM interaction. 



2. The principal axis is odd, and there exists an exact 
real solution with the rotation of I in the easy plane 
determined by the anisotropy. Such an example is 
a system of the type 42^-'2x ^ in which the instan- 
ton trajectory is plane and the DM constant does 
not appear in the real part of the Euclidean action. 
In this case, however, the inclusion of the DM in- 
teraction leads to the emergence of the imaginary 
part of the Euclidean action and may effect on the 
tunneling probability due to the interference of in- 
stanton trajectories lying in different planes. Since 
the corresponding phase factor contains the large 
value N, tunneling can be suppressed completely 
due to the destructive interference even for small 
values of the DM constant. 



3. The principal axis is odd, and the simple solution 

(j) = 27rfc/n does not exist. In this case, the vector 
I has both real and imaginary components, but all 
types of DM interaction change only the real part 
of the Euclidean action, this change being small 
in view of the smallness of the parameter d^ /(J (3). 
The imaginary part of the Euclidean action is equal 
to zero and destructive interference effects are ab- 
sent. 



Thus, the only important effect produced by the DM 

interaction is associated with the possibility of the inter- 
ference of instanton trajectories in the case when their 
number is greater than two (an antiferromagnctic parti- 
cle with the easy magnetization axis approximately cor- 
responds to n > 2). This effect can be observed for anti- 
ferromagnets with an odd principal axis in the case when 
the rotation of I on the instanton trajectory also occurs 
through the odd axis. It is associated with the interfer- 
ence of pairs of instanton trajectories lying in different 
planes. Since I is real in this case and all instanton tra- 
jectories are plane, an exact analysis can easily be carried 
out and the description of tunneling is reduced to the ge- 
ometrical analysis described in the section III. 
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